Relativistic Structure of the Deuteron: 
1. Electro-disintegration and y-Scaling 



(N 
> 

o 
o 

On 
O 
00 
OS 



o 

=5 



X 



C. Ciofi degli Atti, D. Faralli, A.Yu. Umnikov and L.P. Kaptari * 
Department of Physics, University of Perugia, and INFN, Sezione di Perugia, via A . Pascoli, 

Perugia, 1-06100, Italy. 



C$ • (February 9, 2008) 



in 

(N 



Typeset using REVTeX 



*0n leave from Bogoliubov Lab. Theor. Phys., JINR, Dubna, Russia. 



1 



Abstract 



Realistic solutions of the spinor-spinor Bethe-Salpeter equation for the 
deuteron with realistic interaction kernel including the exchange of it, a, uj, 
p, 7] and 5 mesons, are used to systematically investigate relativistic effects 
in inclusive quasi-elastic electron-deuteron scattering within the relativistic 
impulse approximation. Relativistic y-scaling is considered by generalising 
the non relativistic scaling function to the relativistic case, and it is shown 
that y-scaling does occur in the usual relativistic scaling variable resulting 
from the energy conservation in the instant form of dynamics. The present 
approach of y-scaling is fully covariant, with the deuteron being described by 
eight components, viz. the 3 S' 1 ++ , 3 S", 3 D+ + , 3 D^~, 3 P+~ 3 Pf + , 1 P^~, 
1 P{ + waves. It is demonstrated that if the negative relative energy states 
1 P±, 3 Pi are disregarded, the concept of covariant momentum distributions 
N(po,p)j with po = Mo/2 — ^p 2 + m 2 , can be introduced, and that calcu- 
lations of electro-disintegration cross section in terms of these distributions 
agree within few percents with the exact calculations which include the 1 P±, 
3 Pi states, provided the nucleon three momentum |p| < IGeV/c; in this 
momentum range, the asymptotic relativistic scaling function is shown to 
coincide with the longitudinal covariant momentum distribution. 

I. INTRODUCTION 

The concept of scaling, which plays an important role in the investigation of hadronic 
structure, can be introduced in the description of scattering processes whenever the cross 
section factorizes into a product of two different quantities, with the first one reflecting the 
nature of the scattering process, and therefore depending upon the relevant independent 
kinematical variables, and the second one (the scaling function) reflecting the internal struc- 
ture of the target, and therefore depending upon a new variable (the scaling variable) which 
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can be associated to the dynamics of the constituents of the target. 

The best example of scaling (x-scaling) is provided by inclusive Deep Inelastic Scattering 
(DIS) of leptons off nucleons in the Bjorken limit [y — > oo, Q 2 — > oo, xsj = Q 2 /2mu = 
const, with v and Q 2 being, respectively, the energy and the four- momentum transfers, and m 
the nucleon mass), the quantity F 2 (xBj, Q 2 ) = vW^fasj, Q 2 ) where W^(xsj, Q 2 ) represents 
the deviation of the inclusive cross section from scattering off a point-like nucleon, becomes 
Q 2 -independent, i.e. scales in the variable xsj (the Bjorken scaling variable) which can be 
associated to the momentum fraction of the quark inside the hadron. Inclusive scattering of 
leptons off a nucleus A in the Quasi-Elastic (QE) region [y < Q 2 /2m, i.e. XBj > 1) has been 
theoretically shown to exhibit another kind of scaling, the so called y-scaling, which 

can be summarized as follows^: at sufficiently high values of the three-momentum transfer q, 
the quantity |q|PE{| 2 -)(i/, q 2 )/iy^ 2 q 2 ), where the nuclear structure function W{fa{v, q 2 ) 
represents the deviation of the cross section from scattering off a point-like nucleus, scales 
to a function of the variable y, according to \q\Wi^2)(^,q^)/W^(i , ,q 2 ) — > F(y) where, 
in the case of the deuteron (but not for complex nuclei ||), F(y) represents the nucleon 
longitudinal momentum distribution f(y) Q : 

oo oo 

f{y) = J n(p||,p±)dp± = 27r J n(p) |p|d|p|, p\\ = \y\- (1.1) 
Ivl 
Experimental data, due to the effect of the final state interaction (FSI), exhibit only a 

qualitative scaling behaviour, and a quantitative analysis || of the deuteron data [§] taking 

into account (FSI), allowed one to obtain the nucleon momentum distribution in the deuteron 

using eq. (|1.1|) . It should however be stressed that the approach of ref. 0, leading to eq. 

( |1.1|) , is based on a fully relativistic instant-form treatment of kinematics, but on the usual 

non relativistic (Schrodinger) treatment of the deuteron, and therefore lacks of a consistent 

1 For exhaustive reviews of y-scaling see refs. 0, ||, || an d ||. 

2 We will consider, from now on, negative values of y for which the effects of non-nucleonic degrees 
of freedom are kinematically suppressed 
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covariant treatment of the process []. Recently |L1], a relativistic covariant model, based 



on the light-front dynamics and light cone kinematics ||12|| , has been adopted to analyze 
inclusive QE scattering off the deuteron, treating the latter as a system of two spinless 
particles interacting via a simple scalar interaction. Within such a model, the deuteron 
wave function has only one component (the S- wave), the square of which defines the model 
momentum distribution. However, it is well known that in the realistic case of two interacting 
spinor particles, the deuteron state is determined by at least four components if one nucleon 
is on mass shell ||13|| , five components within the spinor light cone formalism |L4], or even 
eight components within the covariant Bethe-Salpeter (BS) approach fl5| , |l6fl . 

The necessity of more than two components (5* and D-waves) in the description of the 
deuteron in a covariant approach, follows from an accurate account of the contribution of 
the negative relative energy states in the deuteron, the so-called P-waves. Frequently, in 
the calculation of observables within covariant formalisms with realistic interactions, the 
final results are rather cumbersome (see, for instance, refs. [|T~3| , p~7| ? p~8|| ) and a separation, in a 
compact form, of an analog of the momentum distribution of the deuteron becomes difficult, 
if not, sometimes, impossible. The inclusive quasi elastic cross section and the concept of y- 
scaling have not been so far considered within covariant approaches with realistic interaction, 
and an investigation of the possibility to define relativistic scaling functions and momentum 
distributions is still lacking. 

In this paper we focus on a detailed study of quasi elastic eD scattering, and the pos- 



3 Note that FSI have been treated in ref. || within a full Schrodinger equation approach, i.e. using 
ground and continuum eigenfunctions of the same Hamiltonian. Such an approach is a correct one 
at the kinematics of the data of Ref. Q which, at negative values of y, correspond to a very small 
(less than the pion threshold) relative energy of the neutron-proton pair in the continuum; it should 
always be kept in mind, however, that at high values of Q 2 , such that the nucleon-nucleon (NN) 
cross section becomes strongly absorptive, the Schroodinger approach is inadequate (see e.g. fllOf ) 
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sibility to analyze it in terms of scaling functions and momentum distributions, using the 
numerical BS solution recently obtained with a realistic one-boson exchange (OBE) inter- 



action 19-21 



Through this paper, as in fllTfl , relativistic effects will be investigated within the Impulse 
Approximation (IA), which means that the final state interaction of the np pair in the 
continuum will be disregarded, though it has been shown || that FSI lead to sizeable scaling 



violation effects at low values of |q|. FSI effects will be analyzed in a subsequent paper |22 
here we are only interested in the investigation of ^/-scaling within a fully covariant treatment 
of inclusive eD-scattering in IA, within a description of the deuteron in terms of realistic 
solutions of the BS equation. 

Our paper is organized as follows: in Section |IJ the derivation, within the BS formalism, 
of the two basic quantities which are necessary to define relativistic ^/-scaling, i.e. the cross 
section for elastic electron scattering from a moving and off-mass-shell nucleon, and the cross 



section for inclusive quasi-elastic scattering from the deuteron, are presented; in Section [D I 
the relativistic scaling function is defined, its non relativistic reduction is illustrated, and 
the results of numerical calculations are presented; a relativistic momentum distribution 
appropriate to the BS approach is defined in Section [TV], where the relationship between the 
relativistic momentum distribution and the scaling function is illustrated; the Summary and 
Conclusions are presented in Section [V|; some relevant details concerning the construction 
of the Mandelstam vertex for the operator of eD scattering and for the computation of 
matrix elements within the Bethe-Salpeter formalism are given in Appendices [A] and ||, 
respectively. 

II. THE ELECTRON-NUCLEON AND ELECTRON -DEUTERON CROSS 
SECTIONS WITHIN THE BETHE-SALPETER FORMALISM 



A. General formulae for the cross section 



In this Section the cross sections for the electron-hadron scattering within the covariant 
BS formalism will be derived. In particular, the elastic scattering from a moving and off- 
mass-shell nucleon, and the inclusive quasi-elastic scattering off the deuteron at rest, will be 
considered. Both processes will be denoted A(e, e')X, where A stands for the target hadron 
and X for the final hadronic states. The 4-momenta of the initial and final electrons in 
the laboratory system are k = (£, k) and k! = {£' , k'), respectively; the four momentum 
transfer is q — k — k' — (v, q), and the orientation of the coordinate system is defined by 
q = (0, 0, q z ). At high energies the electron mass can be disregarded, so that 

k 2 = (k'f^0, kk' = -kq = ^f = ^, (2.1) 
g 2 = - g 2 = 4^ , sin 2 ^. (2.2) 

where 9 is the scattering angle. The following relations will be useful in what follows: 



2 V S1I1 2 / 

cos6» fc = — F J= ( 1 + , (2.4) 



|q| = k 2 | = V^ 2 + ^' ( 2 - 5 ) 

where 6k is the polar angle of the initial electron. 

Hereafter the electron-nucleon vertex in the on-mass-shell form, viz: 

rf (Q 2 ) = 7^i(Q 2 ) + *^(Q 2 ), (2.6) 

will be used, where F 12 are the electromagnetic form factors of the nucleon, and k its 
anomalous magnetic moment. It is well known that the choice ( |2.6j ) for T e ^ N (Q 2 ) violates 
gauge invariance; this is a relevant point which will be briefly discussed in Section |TJ. In 
the one-photon exchange approximation the general formula for the invariant cross section 
for the process A(e, e')X has the following form: 



da 



\M 



4kp A '--- e +^ e " v 1 ™ - ™> (27r) 3 2f 

where is the initial target momentum, the total momentum of the final hadron (X) 
state, drx the phase factor corresponding to the final hadron state X, and A4 e+ A^ e '+x 
the invariant matrix element describing the process. For the elastic electron nucleon (eN) 
(Fig.^a) and electron deuteron (eD) scattering in the impulse approximation (Fig.|T|b), we 
have, respectively: 



, +x \ 2 (2n)^(k + p A -k'-p x ) 



dk' 



-dr 



x ■ 



(2.7) 



Pa = Pi, Px = p'i =Pi + q, 
dp'i 



dr 



x 



(2n) 3 2E[ 



> for eN scattering. 



(2.8) 



p A = Pd, Px = p[ + P2 = Pd + q, 
dp[ dp 2 



dr 



x 



2iry2E[ (2nf2E 2 



> for eD scattering. 



(2.9) 



Using the identity: 



dpi 



- = * ((Pi 



.2 2\ a Pi 



— m 



(2.10) 



(2ti) 3 2E[ ' (2tt) 3 ' 

the elastic eN and the inclusive eD cross sections are obtained by integrating over the d 4 p[: 



da 



\M e+A ^+x\ 2 (2tt)5((pi + qf - m 2 )- 



-dr. 



(2.11) 



dE'dSlw ^kp A J — ' v ' ■*> - / 2(2tt) 3 

where dr is the phase factor corresponding to the final hadron state X without the hit 
nucleon, is the scattered electron solid angle and pi the initial nucleon momentum in 
eN {p\ = m 2 ) and eD (pi = P D — p 2 , pi 7^ m 2 ) scattering. 

The square of the invariant matrix element in eq. ( |2.11| ), averaged over the spins of the 
colliding particles and summed over the spins of the scattered particles, can be cast in the 
form: 



\M 



e+A-^e'+X I 



|2 fn _\ r//_ , _\2 2\ e 



(2n)5(( Pl + qf - m 1 



Q 



■ir(k,q)W£(p A ,q) 



(2.12) 



where the leptonic tensor has the familiar form: 



L" u (k, q) = 2 - (k^q u + k„q„) + j . (2.13) 

As for the hadronic tensor W^ u , appearing in elastic eN scattering, one has 
W^(q 2 , Pl ■ q) = ±Tr {(p, + m)Yf{Q 2 ){p 1 + q + m)Tf (Q 2 )} (2tt)*((pi + q) 2 - m 2 ) . (2.14) 



By contracting (|2.14|) with the leptonic tensor, eq. ( [2.13|) , one obtains, in the nucleon rest 



system ( pi = (m, 0)), the well known Rosenbluth cross section. In the case of a complex 
system, e.g., the deuteron , only the general expression for the hadronic tensor can be 
unambiguously defined in terms of the two independent structure functions, W\^{q 2 ,Pd • <?), 
whose explicit form, however, relies on particular theoretical models; the model adopted in 
this paper is described in the next subsection. 

B. The hadronic tensor 

Our strategy in computing the hadronic tensor for a composite system is the following 
one: a nucleonic tensor operator 0^ u will be defined, whose expectation value between 
relativistic hadronic states | A ) generates the corresponding hadronic tensor, according to: 

W£,= (A\6Z,\A). (2.15) 

The general requirements for the operator 0^ v are as follows: 



i) it should lead to eq. ( |2.14|) when sandwiched between free nucleon states (\A}= \ N)); 

ii) when sandwiched between deuteron states ( | A) = | D)) it should incorporate the effects 
from the Fermi motion and the off mass shellness of the hit nucleon; 

The operator 0^ u , eq. ( |2.15| ), due to the choice of the vertex (|2.6|), has the following 
form: 

6^( Pl ,g) = (27r)5((p 1 + g) 2 -m 2 ) x 



F 2 AQ 2 )6^q) + ^Fi(Q 2 )F 2 (Q 2 )0^( Pl ,q) + —F*{<?)6$fa, q) , (2-16) 



where 



0$(Pi, Q) = 2 [-9»v(Pi + q~m)+ 7^ + 7„p lA4 ] , 
0& 2) (pi, q) = 4 + (g 2 + gpi))] , 

0$(pi, q) = 2 [v(m? 2 + Pig 2 - <?(g 2 + 2gpi)) - q 2 (l^Piu + luPi^) 



(2.17) 
(2.18) 
(2.19) 



and all terms proportional to g^j,) have been omitted in view of the gauge invariance of the 
leptonic tensor ( |2.13| ), q^ u ) LP" = In eqs. fl2.17p -( PTT3D and in the rest of the paper the 
short-hand notation p will be used for the scalar product of a four-vector p with 7-matrices, 



i.e., p = 7 At p At . In actual calculations we first contract Ojj„ with the leptonic tensor L^, and 
the resulting operator is then sandwiched between target ground states. The result of the 
contraction is (see Appendix 0): 



6( Pl , q, k) = L» v 6^ = {6 stat + 50 mot + 50 off } (2tt)5(( Pi + q) 



2 2 
— m 



(2.20) 



where 



O s tat — 2 
6d mot = 4 

so off = 



q 2 m - Amu£ + Am£ 2 A{Q 2 ) + —B{Q 



ni 



mvE — q(kpi) — £{2m£ — rav) + k(2kp\ — qp\) A(Q 



-2q 2 [q+f)F?(Q 2 ) + 2 -^- 
\ 2m / m 



2J1 



(2.21) 
(2.22) 
(2.23) 



q(q 2 + 2q Pl ) F 2 (Q 2 " 



and 



a(q 2 



^ 2 (Q 2 ) - 



-2 9 



^(Q 2 ) 



£(Q 2 ) = [F 1 (Q 2 )+kF 2 (Q 2 



(2.24) 
(2.25) 



Let us discuss the meaning of the terms (|2.21|) -( [2.23|) : O stat represents the contribution 
from the interaction of the lepton with the nucleon at rest, and in case of eN scattering its 
average over nucleon spinors yields the Rosenbluth cross section; 50 mot originates from the 
motion of the nucleon; SO Q ff takes into account the off mass shellness (pf 7^ m 2 ) of a bound 
nucleon, and therefore only apperas in the e-D-cross section. 



As already pointed out, the eN and eD cross sections are obtained by sandwiching the 
relativistic operator 0(pi,q,k) between relativistic nucleon or deuteron states; this means 
that, within the present approach, the hadronic tensor has exactly the same form for a free 
or for a bound nucleon. The only assumption we make, is that the electromagnetic vertices 
for free and bound nucleons are of the same form, with all nuclear effects provided by the 
state vectors. In the usual, non covariant PWIA, based on the use of non relativistic wave 
functions, the cross section off a nucleus A is obtained by relating the nuclear hadronic 
tensor W£ v to the nucleon hadronic tensor W^ u by a convolution formula, with resulting 
ambiguities as far as the extrapolation of the eN cross section for a bound off-mass shell 
nucleon is concerned (see e.g. p3|j24fl ). 

As previously mentioned, gauge invariance is broken for the deuteron current correspond- 
ing to the choice Eq. ( 2.16|) . Several phenomenological prescriptions have been suggested 



to restore it p3|p^| . It should be pointed out that our paper is mainly aimed at theo- 
retically comparing relativistic and non relativistic deuteron momentum distributions and 
^/-scaling functions, without presenting any comparison with experimental data which would, 
of course, require a serious consideration of gauge invariance violation effects. 

The procedure used in this paper can, in principle, be adopted for the description of 
eD scattering within non relativistic Schrodinger picture. However, in this case, consis- 
tency would require a non relativistic reduction of the 'jNN vertex. A systematic study of 
deuteron electro-disintegration within the non relativistic approach taking into account FSI, 



relativistic corrections and meson exchange currents, can be found in |2fJ and references 
therein quoted. 

C. The elastic electron nucleon cross section 

The elastic eN cross section resulting by sandwiching O sta t and 50 mot between free 
nucleon states reads as follows: 

— -^frecaU^, (2-26) 



dVt k > {pxk) 



i n 



where: 

i) f recoil is the recoil term 

! _ 2£sin 2 fl/2-z/ p lz v - 2S sin 2 9/2 

J recoil 

... SEa 



Epi+q I'll Epi+q 



[2.27) 



J P\ 



(pifc) 



comes from the redefinition of the incident flux for a moving nucleon; iii) a eN is 



the "reduced" electron-nucleon cross section, i.e. the cross section without the flux factor 
and the recoil contribution, namely 

,2 



O eN = (TMott 



m 



Epi^pi+q 







A(Q 2 ) + ^ tan 2 -B(Q 2 ) + 



Q 2 

— van -r>^ ^- 2 ££ W /2 



(pik) 2 -( Pl k)( Pl q)-m 2 SS' , (2.28 



where 



0~Mott = 



2 2 f 
CT COS" 1 ; 

AS 2 sin 4 



(2.29) 



is the Mott cross section. In the above equation E, 



pi 



m 2 + pf , J5j 



+ (pi + q) 2 . It can easily be seen from ( |2.26 )-( 2.28 ) that for a nucleon at rest 



(pi = (m, 0)) eq.( [2.26| ) coincides with the Rosenbluth cross section. 



D. The inclusive electron deuteron cross section within the BS formalism 



The relativistic impulse approximation for the inclusive e-D-cross section is obtained by 
averaging the operator O, given by eq. (|2.20| ), with the Bethe-Salpeter amplitudes (see 
Fig. H). The result is 



da 



BS 



Pmax 



° Mott 



pMpI 



2tt 



d(p p 



1 



(2vr) 4 16££' cos 2 f J \p\-\ci\2wM D (M D (M D -2w)) 2 

y\ 



d£'dQ k , J eD 

o Tr {Gm(p, Pd)(pi + m) (6 s tat + SO mot + 50 o //) (pi + m)G M {p, Pd){pi + m)} 



(2.30) 



where p\ = Pd — P2, Gm is the Bethe-Salpeter D — > iVW vertex, M is the deuteron total 
angular momentum projection, p is the relative momentum of the nucleons, i.e. px,2 = 

1 1 



Pd/2 ± p. The deuteron-nucleon vertex Gm is the truncated Green's function which is 
related to the conjugated Bethe-Salpeter amplitude by 

* w(p ' p ^ = Pd) m^) • (2 - 31) 



The limits of integration in ( |2.30|) , |p| m j n = \y\ and p maa ; are obtained from energy 



conservation provided by the 5- function 



5{{P D -p 2 + q) -m ) = 5 [cos 9 P , (2.32) 

2 P ■ q V 2 p • q ' 



which determines the value of cos6 p from the constraint: 



-1<cos0 p <1. (2.33) 



Solving the inequalities ( p.33| ) with cos 8 P defined by the argument of the 5-function ( p.32| ) 
we obtain the same result as in [Pfl, i.e. 



iPl 



mm 2 



/ Am 2 
{M D + u)\l |q| 



\y\ (2-34) 



1 / 4m 2 

|p|max = - I (M D +V)]ll — + |q| } = Pmax, (2.35) 

where s denotes the Mandelstam variable for the 7*-D vertex 

s = (P D + qf = M D (M D + 2u) - Q 2 . (2.36) 



It can be seen from (|2.35|) that p ma x sharply increases with |q|; such a circumstance, as we 



shall see, has relevant consequences for the occurrence of y-scaling. 

In the calculation of the trace appearing in eq. (|2.30|) , two different operators have to 
be considered, namely the scalar operator, 1, coming from O s t a t{ eq. (|2.21|) ), and the vector 
operator 7^, contained in 50 mot and S0 ff, (eqs. ( |2.22|) and ( |2.23| )). 



We have calculated the BS average of these operators, {i) po i e and (j^poie, where the 
subscript "pole" means that the trace in eq. ( |2.30| ) has been obtained using eq. ( |2.31| ) and 



by evaluating the integral over p 2 o in the pole corresponding to the second particle on mass 
shell (for details, see Appendix IB), i.e., 



1 9 



P20 = W 



p 2 + m 2 , p 



M, 



D 



- w, p w = M D -w. 



(2.37) 



The eD-cross section can now be rewritten in terms of the nucleon pole contributions to 
the vector and the scalar parts, namely 



da 



dS'dQk' 



BS 



eD 



Pmax 



\p\ 2 d\p\ 



2tt 



\u\ 



(2tt) 4 8££' cos 2 f { 2 IpI ' 1^1 



dcpp 



(fstat + Sfmot + Sfoff) , (2.38 



where 



2m 



f stat 

Sfmot = 4 
8 faff = ~2q 2 



q 2 + A££' A(Q 2 ) + ^B(Q 2 



m 



(i) P B i(p), 



-2m££' (i>*f e (p) + (-q»(k Pl ) + k»(2k Pl - q Pl )) ( 7 „>Se(p)l MQ 



[2.39) 
(2.40) 



\ BS 

2nq 2 



m 

2„2 



pirn 

mq^)?l(p) + qpi(l&)} F 1 (Q 2 )F 2 (Q 2 ) 



k q 



2m 2 



q"(iXie(p)^ + 2m)}F 2 (Q 2 ). 



(2.41) 



Let us compare equation ( |2.38| ), which represents the BS eD inclusive cross section in 
which the matrix elements of both the eN and D — > iViV vertices are treated covariantly, 
with the usual non covariant PWIA cross section, where the eN vertex is treated covariantly 
within the instant-form of dynamics, and the vertex D — ► NN is treated non covariantly 
within the Schrodinger approach, viz. 



da 



d£'dQk' 



PWIA 



eD 



Pmax 



(27T) 



|p| d\p\n D (p[ 



E, 



p+q 



S^ep @ en) i 



(2.42) 



\y\ 



where a ep ^ is the relativistic electron-nucleon cross section for a free moving nucleon (see, 
e.g., ref. |23 ) and the non relativistic momentum distribution n D (p) is normalized as follows: 



J dpn D (p) = 1. 



(2.43) 



The most relevant difference between eqs. fl2.38|) and ( |2.42| ) arises from the non covariant 
treatment of the D — > NN vertex which, in eq. ( |2.42| ), is entirely determined by a single 
quantity, the nucleon momentum distribution rao(p), whereas in the covariant BS cross 



i 3 



section, eq. Q2.38Q , it depends upon (iV£f e and (jn)^ e - In order to exhibit the quantitative 

)pl and ( 7 „)££ e 



difference between the two cross sections let us compare (i)£f e and (7^)p f e with their non 



relativistic limits (|p| /m « 1) which are given by (for details see refs. [19,26]): 



(7M>3e(p) - W(P,/E P ) n D (p), (l}^ e (p) - (27r)^(m/E p ) n D (p), (2.44) 

where = (£p,p) and no(p) is the deuteron momentum distribution. 

In Figs.|3] and £| (l)™f e and (7 M )p f e are compared with their non relativistic limit ob- 
tained with nz)(p) corresponding to various realistic interactions. Using eqs. ( j2.44 ) the non 



relativistic limit of the BS cross section can be obtained straightforwardly, viz. 

, \ NR Pmax 

(2tt) f \ p \d\p\n D {p)-^f{a ep + a en ), (2.45) 
ii 111 



\dS'dfl k r J eD 

where a e N = or e jv(|p|, |q|, v) is given by eq. ( [2.281 ), without the off mass shell contribution 
(eq. ( |2.41|) ), which can easily be shown to vanish in the non relativistic limit. 

Thus, we have demonstrated that the non relativistic limit of the BS inclusive cross 
section (|2.38|) , obtained by taking the non relativistic limit of (lV£f e and (7 M )p G f e , has exactly 



the same structure as the instant form result of ref. 0], apart from some minor differences 
between the relativistic eN cross sections a e N and a e N, which are irrelevant for the present 



paper, and which will be discussed elsewhere [22 



In closing this Section, the following remarks are in order: 

i) the BS covariant inclusive eD cross section does not factorize into a product of an electron 
nucleon cross section and a deuteron structure function. In that respect the covariant results 
differ from the usual non covariant PWIA; 

ii) within the BS formalism the interacting nucleon is consistently treated as an off mass 
shell particle. Consequently, the matrix element of the 'jNN vertex is half off shell for the 
eD scattering. As a result, additional off mass shell effects, represented by eq. (|2.41|) , arise 
due to covariance of the approach. 



III. THE RELATIVISTIC SCALING FUNCTION 



A. Non relativistic and relativistic scaling functions 



In the non relativistic case, the concept of ^/-scaling can be introduced when the value 
of |q| becomes large enough so has to make p ma x ~ °o an d the dependence of cx eA r ( or a eN ) 
upon |p| very weak. In such a case eq. (|2.45|) can be cast in the following form: 



da \ NR , , E, 



d£'dVt h > 



Sep + S en ) (2tt) J | P | d\p\n D ( P ), (3.1) 

Ivl 



K / eD 

where s e N and -E^+ig] represent a e N and -Ep+g- , calculated at |p| = |p| m j n = \y\ and are taken 
out of the integral. Such an approximation has been carefully investigated in ref. |J and 
found to be valid within few percents, provided Q 2 > 0.5GeV 2 /c 2 . It is clear therefore, that 
at large values of |q| the following quantity (the non relativistic scaling function) 

II / rl \ NR 

will be directly related to the longitudinal momentum distribution || 

F NR (h\,y) — » f(y) = 2tt| IpMpK(IpI), (3-3) 

\y\ 

whose first derivative will provide the non relativistic momentum distribution. As already 
pointed out, the condition for the occurrence of non relativistic ^/-scaling is that eq. (|2.45|) 
could be cast in the form ( |3.1| ), which means that: i) Q 2 > 0.5GeV 2 /c 2 , in order to make 
the replacement a e N — > s e N and E$ + $ Ey+\$\ possible, and ii) p ma x = (|q| — \y\) \y\ 
(cf. (|2.34|) and (|2.34|) ) in order to saturate the integral of the momentum distribution, 

Pmax CO 

j |p|d|p|n£)(|p|) — > j |p|d|p|n£)(|p|). Condition ii) obviously implies that the larger the 

Ivl \y\ 

value of \y\, the larger the value of |q| at which scaling will occur. The satisfaction of 
the inequalities |q| ^> 2\y\, XBj > 1 leads, for any well-behaved TO£>(|p|), to the following 
conditions for the occurrence of non relativistic ^/-scaling: 

2m/3 < v < |q|, |q| > 2m. (3.4) 

Note, that the above conditions are very different from the conditions for Bjorken scaling 
v ~ |q|. Let us now discuss relativistic scaling. To keep contact with non relativistic scaling, 
let us define the following relativistic scaling function: 

1 R 



F BS (\dly) = #- • i-J^A BS I (s ep + s en ) , (3.5) 



eD 



, \ BS 

da 



with given by eq. ( |2.38|) . In the rest of the paper the following questions will 

\ dc d\ l k i J eD 

be dadresse: 

i) does (and at which values of |q|) eq. ( p75| ) scales in y ? 

ii) if scaling does occur, can a relationship be established between the asymptotic scaling 
function and the momentum distribution? 



It is clear, by looking at eq. Q2.38Q , that relativistic scaling, as the non relativistic one, is 
entirely governed by the value of p ma x (eq. (|2.35| )), in that one expects that, starting from 
a certain high value of |q|, p max becomes large enough as to saturate the |p| dependence of 
(l)po/e(|p|) and (7 M ) po ze(|p|). It should be pointed out that if scaling of ( |3.5|) is observed, this 
would imply that a factorisation of ( p. 38 ) similar to the one occurring in the non relativistic 
case (i.e. eq. ( |3.2|) ) occurs in relativistic scaling as well. Such a factorisation, due the 
complex structure of eq. ( j2.38j ) is, a priori, not obvious. As for the second question, it is 
by no means a trivial one, for in the BS case even the concept of momentum distributions 
is not well defined. Nevertheless, we will see that the concept of relativistic momentum 
distribution can be introduced, and that a relationship of such a momentum distribution 
with the asymptotic scaling function, can be established. 



B. Numerical calculations of the relativistic cross section and scaling function 

In this section the results of numerical calculations of the relativistic scaling function 
F (|q|, y), eq. (|3.5|) , will be presented . In our calculations the numerical solution [[i9ipl| of 



the spinor-spinor BS equation containing a realistic one-boson-exchange interaction kernel, 
which includes the set of tt, a, u, p, rj and 5 exchanged mesons, is used. The meson 
parameters (masses, coupling constants and cut-off parameters) have been taken to be the 
same as in ref. |15.[16|], except for the coupling constant of the scalar a-meson, which has 



been adjusted to provide a numerical solution of the homogeneous BS equation. Recently, 



the solution became available in the form of analytical parametrizations obtained by fitting 
the numerical solution to the BS equation, using the least-squares procedure pOp . 

The details of the numerical calculation of various matrix elements appearing in eD 
electro-disintegration are given in Appendix [FJ. 

In Fig.|5| the approach to scaling of the BS scaling function ( |3.5| ) is shown for various 
values of y. It can be seen that scaling is approached very rapidly due to the sharp increase 
of Pmax in eq. ( j2.38 ) with |q| (cf. Fig.|6]). It can also be seen that the value of |q|, at 



which scaling is reached, sharply increases with the value of y, going from |q| ~ lGeV/c 
( v ~ 0.3GeV ), at y = -0.2GeV/c, to |q| ~ 2GeV/c {y ~ 0.8GeV) at y = -0.8GeV/c. 
These values match very well the condition for non relativistic ^-scaling (|3.4j) ; this apparently 
surprising result will be explained later on. Let us now briefly discuss the role played by 
inelastic channels on the scaling function. It can be shown that for all values of y 
presented in Fig.[5|, the values of v and |q| in the region of the approach to scaling are below 
the pion production threshold, which can be reached only at higher values of |q|, provided 
Q 2 < 5 GeV 2 /c 2 . Therefore, inelastic excitations of the nucleon are kinematically forbidden 
in a wide range of |q| for a given value of y. The asymptotic scaling function is shown 
vs. y in Fig.[7], whereas the contribution of the off mass shell terms ( [2.41| ) is presented in 



Fig|| By comparing the two figures, it can be seen that the off mass shell corrections are 
negligibly small in the whole range of \y\, and in what follows they will be disregarded. 
Figure |9] illustrates the role of the "moving" components, eq. ( |2.40|) , calculated at different 
values of |q| = 3,10 and 18 GeV/c; as expected, these corrections increase with \y\ and 
are practically |q|-independent. The various contributions to the total scaling function are 



presented in Pig.|Ty, whereas, in Fig.|TT] the asymptotic relativistic BS scaling function is 
compared with the non relativistic one, calculated with various realistic interactions. It can 
be seen that for \y\ > 0.4 GeV/c, the differences between the BS and the non relativistic 
scaling functions are very large, except for the RSC interaction. 

The scaling behaviour of the relativistic scaling function (|3.5| ) shown in Fig.[5] would 



imply the factorisation of the BS inclusive cross section (|2.38|) into the free eN cross section 
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and some kind of deuteron structure function. Due to the complexity of eqs. (|2.39|) - ([2.41|) , 
neither the origin of such a factorization, nor the nature of the deuteron structure function 
are clear at the moment; they will however be discussed and clarified in the next Section. 



IV. RELATIVISTIC MOMENTUM DISTRIBUTION 

Since in the covariant BS formalism the deuteron amplitude does not have a probabilistic 
interpretation, the concept of momentum distribution is ambiguous. Nevertheless, we will 
now rearrange the matrix elements of 1 and 7^ in such a way that, under certain conditions, 
they could be interpreted in terms of a relativistic momentum distribution. 

Let us return to the main quantity, eq. fl2.30| ), and try to analyze it analytically in 
more details. To this end, it is convenient in the decomposition of the BS amplitude, to 
shift from the Dirac basis, used in ref. |L9], to a basis of spin-angular matrices [|15| , |30|] , 



i.e. an outer product of two spinors, representing the solutions of the free Dirac equation 
with positive and negative energies. This basis, which is frequently used, is labeled by the 
relative momentum p, the helicities Aj and the energy spin p, of the particles [|L5|| , and is 



sometimes called the (J, Ai, A2, Pi, P2) representation. The spectroscopic notations are used 
for the partial amplitudes viz. 2S + 1 £ J p J 1 ' p2 ^ j. e-) 

3 S++, 3 S r , 3 D+\ 3 D r , "Pf, ^f + > 3 Pi + - (4-1) 

The partial amplitudes in the basis ( [4.1|) exhibit a more transparent physical meaning, since 
they can be compared with the deuteron states in the non relativistic limit. It is intuitively 
clear that the two nucleons in the deuteron are mainly in positive energy states with L = 0,2, 
so that one may expect the probability of negative energy states with L = 1 in eq. (|4.1|) to be 
much smaller than the probability for the 3 S^ + and 3 Df + states. Moreover, it can be shown 
that the waves 3 Si + and 3 Df + directly correspond to the S and D waves in the deuteron, 
with the waves with negative energy vanishing in the non relativistic limit. (The connection 
between the partial amplitudes defined in the Dirac and the p-spin representations can be 
found in ref. jTBfl). 
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Let us now investigate analytically the matrix elements ( 7 M )£f e (p) and ( 1 )p Q f e (p), eqs. 
( 2.39 )-( 2.41j) . Their explicit expressions are: 



— z- 



6M 



D M 



(pi — m 2 ) 2 (p 2 — m 2 ) 



Tr 



G M (pi + m)6 (pi + m)G M (p 2 + m) 



2tt 



(4.2) 



pole 



where O stands for either 7^ or 1. Now instead of calculating the pole contributions in 
eq. ( |4.2j ), we go back to the original definition of these averages, following eqs.( [2.30D -( [2.31D , 
namely 

/ (p 2 — m 2 ) 1 



^o(p) 



J^Tr ^ M 0^ M {p 2 -m) 



(4.3) 



p =M D /2-E 2 



\ AM D w 3 hi 

and calculate directly the trace (|4.3|) evaluated at p = M D /2 — E 2 . Here it is worth 
emphasizing that in our case, when one particle (the second one in the present notation) is 
on mass shell, only four partial amplitudes contribute to the process [O] , namely, only those 
partial amplitudes (|4.1|) with the second p-spin index positive, i.e. the 3 Si + , 3 Df + , 1 P{ + 
and 3 P 1 _+ amplitudes. In correspondence to these contributions, we introduce the Bethe- 
Salpeter wave functions for each vertex, viz. 



^s(po, IpI) = 
* ft (po,|p|) 



1 G+ + ( P o,|p|)/27T 

V2n V2M^(2E 2 - M D ) ; 

1 gtZSpoM)/**. 



^d(po, IpI) ^ 

*ft(Po,|p|) 



1 



g5 + (P0,|p|)/27T 

V2tt v / 2Md(2E 2 - M D ) 
1 G 3 -+(po,|p|)/2vr 



(4.4) 



(4.5) 



where the normalization factors have been chosen so as to correspond to the non relativistic 
normalization of the deuteron wave function: 



J dp (u 2 (p)+w 2 (p] 



1. 



Then for TZq we obtain: 



n 10 = (2vr) 3 (*|(po, |p|) + ^(p , | P |) + tt^fo, | P |) + *| 3 ( Po , | P |)) ; 
^ = (2vr) 3 1- (*|(po, |p|) + *f,(po, |p|) - ^(p , | P |) - *^(p , | P |)) + 57^; 



E 2 
rn 
E 2 



Hi = (2tt) 3 ^ (v&|(p , |p|) + * 2 D ( Po , |p|) - ^> , IpI) - 4^(pa, IpI)) + sn h 



(4.6) 

(4.7) 
(4.8) 
(4.9) 



1 



where 



67U = (2tt) 



P [2^3 



rn 



E 2 I 3|p| 



SR 1 = -{2k) s ^- 



in | 2v3|p 
i?2 I 3m 



^s(po, |p|) (*^(po, IpI) - v^*j%(po, IpI) 
*2,(po, IpI) (V2* Fl (p , IpI) + *p 3 (po, IpI))] } 
#s(p , IpI) (*Pi(Pd, IpI) - v^* ft (po, IpI) 



(4.10) 



+*u(po, IpI) (V2* Pl (p , IpI) + ^p 3 (po, IpI))] } ; (4ii) 



As a conclusion, the inclusive eD cross section will be given by eq. (ggp with eqs. 
(fOgD-(CT) for (I) 55 and replaced by eqs. (|47|)- (pl]) . 

From what we have exhibited, it can be seen that in the BS formalism, there is no 
universal momentum distribution (cf eqs. (|4.7|) - ( |4.9|) ) so that, in principle, a factorized 
cross section in the form ( 2.45Q does not hold. This is a consequence of the covariance of the 
BS formalism, where the small components (|4.5|) with negative relative energies are taken 



into account. It has been shown []TSJ that the contribution from the waves with positive 
relative energies are much larger than the one from typ(p , |p|), which therefore, can be 
disregarded. However, the corrections ([4.10 ) and ( |4.11| ), resulting from the interference 
between large and small waves contribute both to the static and to the moving nucleon 
contribution to the eD cross section and, "a priori", cannot be disregarded; accordingly 
they will be taken into account in our calculations. Let us now introduce the following 
quantity, which will be called hereafter the covariant relativistic momentum distribution : 



N S5 (p ,p) = N(p ,p)+5N(p ,p), 



(4.12) 



where 



N(po, p) = (*s(po, IpI) + ^Upo, Ip 

5N( Po ,p) = {^ 



(4.13) 



^s(po, IpI) (*ft(Po, IpI) - V2^ Ps ( Po , | P 
q D ( Po , |p|) (V2^ Pl ( Po , |p|) + ^ P3 {p , IpI) 



(4.14) 



on 



Since the relative energy p is fixed, the momentum distribution ( |4.13| ), which is defined 
only in terms of the S and D components, resembles the non relativistic distribution no(p), 
and therefore it is expected to provide the main contribution to the eD cross section. 
Thus the matrix elements TZq, eqs. ( |4.7| )-( |4~9|) can be written in the following way: 



TZ^ = (2tt) 3 ^ ■ { m (4.15) 

E * N(po,p) + ^5N(p ,p), /i=(l,2,3) 

^ IpI 

Ki = (2vr) 3 ^ ■ {n( Po ,p) - M 5N(p0)P )| . (4.16) 
E 2 { m J 

If, for the time being, the contribution of 6N(po, p) is disregarded , it is possible to relate 
the BS inclusive eD cross section to the elastic eN cross section for a moving nucleon; as a 
matter of fact, by inserting eqs. (|4.15| )- (|4.16| ) into eqs. (|2.39| )- (|2.41|) , one obtains: 



N(po,p), ^ = 0; 



\ BS Pmax , 

da \ , „ r , , ., , hi, 



(2tt) / \p\d\p\N( Po ,p)^f(a ep + a en ), (4.17) 



dS'dQu I „ J |q 

/ eD |y| ! 



At large values of |q| eq. ( |4.17| ) becomes 



da \ r , hi\ 



dE'dVty 



{Sep + S en } ^(2tt) / | P |d| P |N(po, P ). (4.18) 

\v\ 



where, as before, a eN and Ep+Q have been evaluated at |p| = |p m j n | = \y\. If eq. ( |4.18| ) is 
placed in ( p.5|) , the BS asymptotic scaling function is obtained, viz. 

I, / , \BS oo 

f BS (y) = -^-(-^r) ■{s ep + Sen }- 1 = (2n) [ |p|rf|p|N(p ,p). (4.19) 

Ey+\$\ \ dt dil k'JeD /, 

from which information on the covariant nucleon momentum distribution N(p , p) could be 
obtained. 

The asymptotic scaling function, calculated by eq. ( |4.19| ), coincides with the exact 
scaling function (obtained using eqs. (|2.38|) -( [2.41|) in eq. ( p.5|) ) in the whole range of y (the 



largest difference occurring at y — —0.8 and being less than 10%), which means that the 
BS inclusive cross section could be safely replaced by its approximation ( 4.18| ). In order to 



understand such a result, in Figs. |T2|- [TJ] the separate BS waves, viz. the S and D waves 
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(full lines of Figs. and ]13|), and the P\ and P3 waves (Fig. [14] ), are presented. It can 
be seen that the latter waves are smaller by order of magnitudes than the former ones, so 
that the quantity 5N(po,p), due to the interference terms (|4.10| ) and ( |4.11| ) generated by 



the negative energy states, turns out to be negligible up to |p| ~ lGeV/c. Thus we can 
conclude that in the interval < |p| < lGeV/c the total distribution ( |4.12| ) can be safely 
approximated by the diagonal contribution ( |4.13| ), as illustrated by Fig. 15, where the full 



momentum distributions are shown. This is the reason why the BS inclusive cross section 
factorizes in the same way as the non relativistic one, and the relativistic scaling function 



75]) scales in y. In Figs. [L2|- [15] we have also compared the full BS results with the results 
from other approaches, viz. the non relativistic Schrodinger approach with different types 
of interactions (see also ref. @), as well as the relativistic approach based upon the Gross 
equations |L3|. The BS and Gross approaches differ both in the form of the relativistic 
equations as well as in the number of exchanged bosons considered in the kernel (six in the 
former approach and four in the latter one), but both reproduce equally well the experimental 
NN phase shifts and the ground state properties of the deuteron, which is reflected in the 
very similar behaviour of the S and D waves shown in Figs.[12| and these results also show 
that the high momentum content (|p| ~ 0.5GeV/c) generated by the relativistic equations 
is appreciably higher than the one provided by non relativistic wave functions. 

V. SUMMARY AND CONCLUSIONS 

In the present paper the inclusive quasi-elastic electron-deuteron cross section has been 



analyzed within the relativistic Impulse Approximation using recent, realistic solutions |2T 
of the spinor-spinor Bethe-Salpeter equation for the deuteron, with the interaction kernel 
including the exchange of n, a, u, p, r] and S mesons. In our approach, both the 7*iV and 
the D — > iViV vertices are treated relativistically, with eight components for the deuteron 
amplitude, unlike the usual, non relativistic approach ||, in which the j*N is described 
by a relativistic free electron-nucleon cross section, and the D — > iViV vertex by the usual 
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non relativistic two-component Schrodinger wave function (we reiterate that when we call 
the latter approach "non relativistic", we refer only to the D — > NN vertex) . The aim of 
our paper was twofold, viz.: i) to investigate the relevance of relativistic effects, and ii) to 
understand whether the concept of y-scaling can be introduced in a relativistic description 
of inclusive e-D-scattering. The main results of our analysis can be summarized as follows: 

1. The relativistic inclusive eD cross section has been obtained in terms of the pole 
matrix elements (l)p f e (|p|) and (7^)p f e (|p|), taking into account the off-mass shellness of 
the nucleon and it has been found that unlike the non relativistic case, the BS cross section 
does not factorize into a product of the free electron-nucleon cross section and a structure 
factor depending upon the deuteron momentum distribution. 

2. It has been shown that the BS cross section can be written as a function of the three 
momentum transfer q and a variable y, which is exactly the same relativistic scaling variable 
used in the non relativistic approach and resulting from the relativistic instant-form energy 
momentum conservation. Thus, in full analogy with the non relativistic case, a relativistic 
scaling function F BS (\q\,y) has been defined as the ratio of the BS eD cross section to the 
free eN cross section (times a proper phase space factor), and y-scaling of F BS (\q\,y) has 
been demonstrated to occur, i.e. F BS (\q\,y) — > f BS (y), with the conditions for relativistic 
y-scaling being very similar to those of non relativistic scaling, i.e. 2m/3 < v < |q|, |q| > 
2m. 

3. It has been pointed out, that whereas the mechanism of non relativistic scaling is easily 
understood in terms of the rapid decay of the momentum distribution np(p), which makes 

the non relativistic scaling function F NK (\q\,y) ~ 2tt / n(p) |p| d|p| to rapidly saturate 

\y\ 

with |q|, i.e. to scale in y, a similar explanation in the relativistic case is not, in principle, 
possible, since, as stated in point 1, the BS cross section does not factorize, and, moreover, 
the concept of momentum distribution in the BS case is not uniquely defined. Thus, in order 
to understand the mechanism of the observed relativistic y-scaling, the role of the various 
components of the BS amplitude was analyzed, and it has been found that if the extremely 
small diagonal contribution of the negative energy P-waves is omitted, it is possible to define 

9^ 



a covariant momentum distribution of the form N ss (p ,p) = N(po 5 p) + ^N(p ,p), where 
5N(po,p)> which originates from the interference between the positive and negative waves, 
can be safely disregarded provided |p| < lGeV/c, so that, as a result, the BS cross section 

factorizes in the same way as does the non relativistic cross section and the relativistic 

\i\-\v\ 

scaling function becomes F (|q|, y) ~ 2tt j N(p , p) |p| ^|p|; such a result provides the 

Ivl 

explanation for the relativistic y-scaling and makes it possible to obtain the BS covariant 
momentum distributions by a simple first order derivative of the asymptotic BS scaling 

oo 

function f BS (y) ~ 2tt J N(p ,p) |p| d|p|- 

Ivl 

4. The BS relativistic momentum distribution, N(p ,p), and the non relativistic one , 
ri£)(p), are practically the same up to |p| ~ 0.4 0.5 GeV/c, where they start to differ by 
an amount which depends upon the two-body interaction producing no(p). 

To sum up, it can be concluded that, if the effects from negative energy P-states can be 
disregarded, which has been demonstrated to be the case when the nucleon momentum in 
the deuteron |p| < lGeV/c, the concept of y-scaling can be introduced in the BS relativistic 
description of inclusive quasi-elastic eD scattering, in the same way as it is in the conven- 
tional non relativistic approach, i.e. by introducing a scaling function which, in the scaling 
regime, is nothing but the nucleon longitudinal momentum distribution; moreover, both in 
the relativistic and non relativistic cases, scaling is shown to occur in the same variable y, 
and at values of v and |q| such that quasi-elastic scattering is the dominant process. 
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APPENDIX A: THE CONSTRUCTION OF THE NUCLEON OPERATOR 



The contraction of the operators ( [2.16| )-( |2.19| ) with the leptonic tensor ( |2.13| ) is given by: 



d( Pu q,k) = d, u ( Pl ,q)L^(k,q) 



(Al) 



F?(Q 2 )6W( Pl , q, k) + ^F l (Q 2 )F 2 (Q 2 )6^(p 1 , q, k) + £^F 2 (Q 2 )6^( P i, q, k) (A2) 



x(2n)5(( Pl + q) 2 -m 2 ) 



(A3) 



where 



6 (1) (pi, q,k) = 2 q 2 m - q(q 2 + 2k Pl ) + k(Ak Pl - 2q Pl ) , 
6 {12) ( Pl , q, k) = Aq 2 [(-mq + (q 2 + q Pl ))] , 

0^ 2 \ P i, q, k) = 2q 2 mq 2 — q(q 2 + 2q P i — 2k P \) — 2k(2k P \ — qpi) 



(A4) 
(A5) 
(A6) 



The average of (|A6|) with the nucleon or the deuteron amplitudes, times the factor e 2 /Q 4 , 
gives the invariant matrix elements. 



A graphical representation of the operator 0( P \,q,k) is presented in Fig. 16|, where 
the crossed nucleon line corresponds to a nucleon on the mass-shell with the propagator 
(2ir)8(( P i + q) 2 - m 2 ) (pi + q + m). 

The nucleon operator 0(pi, q, k) is the central object in our discussion about the connec- 
tion between the nucleon and the deuteron cross sections. Let us first consider the differential 
elastic cross section for the process e + N — > e' + N. The invariant matrix element is defined 
by 



el 

\M e +N^e'+N\ 2 (2tt)5 (( Pl + qf - 771 2 ) = — - £](pi, S l\0( P l, q, Si) 

e 4 1 r - 

Let us first obtain the cross section in the rest frame of the nucleon ( P i = (m, 0)) 

-Tr {(pi + m)q} = 2mis, -Tr |(pi + m)&j = 2m£, — Tr {(pi + m)} = 2m, 

fcpi = m£, = mi/, 



(A7) 
(A8) 



(A9) 
(A10) 



9K 



so that 



\M e+N ^ +N \ 2 (27r)<y((p 1 + qf - m 2 ) = 

4 (m 2 q 2 - mvq 2 - Am 2 uS + Am 2 S 2 ) F 2 (Q 2 ) + AKq 4 F 1 (Q 2 )F 2 (Q 2 ) + 



e 

Q 1 



K 2 q 2 (q 2 -4S 2 + AuS)F 2 (Q 2 



{27r)5[2mv - Q' 



(All) 



^-16m 2 ^' 
Q 4 



Q ^2®\ -^21^ K< 1 -2 & Tp frflST? t r\V\ 



cos - 

2 2m 2 



sin 2 ~ F 2 (Q 2 ) - ^ sin 2 ^(Q 2 )^ 



N 



(27r)5(2mz/ - Q s 



(A12) 



/) / i £ -2„2 \ „2 /) / 



cos^ - ( Fi(Q 



(27r)5(2mz/ -g 2 ) . 



(A13) 



Inserting the last expression in eq. ( |2.11| ) we get the Rosenbluth cross section: 
da 



a 2 m 



dS'dQw 2S 2 sin 4 i 



x 



cos 2 t {F 2 (Q 2 ) - ^F 2 (Q 2 ) ) - ^ sin 2 v - ( F l( Q 2 ) + nF 2 (Q 2 ^ 



Am 2 



5L 

2m 2 




2 



CMott 



k2( 1 2 p2u\2 , 

Am 2 



^tan 2 ~U(Q 2 )+kF 2 (Q 2 ) 



5(2mv-Q 2 ) = (A14) 



< Ai5 » 



Let us now consider the eD cross section. To this end, we will consider an arbitrary 
reference frame where the four momentum of the nucleon is p\ = (pio>Pi); moreover the 
nucleon can be off mass shell ( p\ 7^ m 2 ). 

Our strategy now is to explicitly separate that part of the operator describing the free 
nucleon at rest, i.e. the part defining the cross section (|A15 ), from the remaining parts of 
the operator due to the nucleon motion and the off-mass-shell corrections. 

Rearranging the terms in eqs. ( |A1[ )-([A^) according to their contributions, we find 

0( Pl , q, k) = {6 on + 50 off } (27T)<y((p 1 + q) 2 - m 2 ) (A16) 

where O on is the sum of terms contributing to the invariant matrix element for the free 
nucleon, and S0 fj is the sum of terms providing the off-mass-shell corrections. Since we 



OR 



can add to and subtract from both terms in eq. (|A16 ) pieces vanishing for a nucleon on the 
mass-shell, the separation on O on and 50 ff is not unique, but the physical results do not 
depend on the way this is done. In our definition: 



O on = — B(Q 2 
m 



q 2 m - 2q(kpi) + 2k(2kp 1 - qpi) A(Q 



50 



ff = 

* I z i ^ 2 i T^2/ r ^2\ , , . , „ 2 \ 77 rr\2\ k2< 1 2 



\ 2m I m 



, -mq + qp 1 )F 1 (Q 2 )F 2 (Q 2 )- 



2m 2 



(A17) 
(A18) 
q(q 2 + 2q Pl )]F 2 (Q 2 ). 



It can be seen that in eq. ([A 171) A(Q 2 ) and B(Q 2 ) are weighted by two different functions, 
viz. a scalar, and a mixture of scalar and vector terms, respectively. This is a result of the 
relativistic structure of the our formalism, in which, unlike the non relativistic case and the 
relativistic case for a nucleon at rest, the vector charge density does not coincides with the 
probability density. Because of the different structure of the weighting factors of A(Q 2 ) and 
B(Q 2 ) in O on , the average value of the latter calculated for the deuteron, will not factorize 
into a common term for A(Q 2 ) and B(Q 2 ). However, if the difference between the vector 
and the scalar charges is not too large, one could be able to define a common term for both 
A(Q 2 ) and B(Q 2 ), plus proper correction terms. Therefore, we redefine O on in the following 
way: 



on = O stat + 50 



'moti 



O s tat — 2 

50 mot = 4 



q 2 m - 4mis£ + Am£ 2 A(Q 2 ) + —B(Q 2 ) 



m 



mvE — q(kpi) — £(2m£ — rav) + k(2kp\ — qpi) A(Q 



(A19) 
(A20) 
(A21) 



where O s t a t is a scalar operator defining the invariant matrix element for the nucleon at rest, 
whereas SO mo t, which is defined by the combination of scalar and vector currents, gives a 
non-zero contribution for the moving nucleon. 
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APPENDIX B: THE POLE STRUCTURE OF THE MATRIX ELEMENTS 

WITHIN THE BS FORMALISM 



In computing the deuteron observables within the BS formalism, i.e. the matrix elements 
of a given operator (D\0\D) = (O), one makes use of the Mandelstam technique, which 
yields 

d A p 1 



(O) 



2M D J (2tt) 4 (pi - m 2 + ie) 2 (p 2 2 - m 2 + ie) 

\j2 Tr {Gm(Po,p)(Pi +m)6(p 1 + m)G M {po, P){p2 + m)} . (Bl) 



3 j j 



It can be seen that in eq. ( plf ) there are poles and cuts in the integration variable pq. 
However, the whole matrix element (O) is real and finite. This allows one to perform the 
Wick rotation in the complex plane of the relative energy p , and to safely compute the 
integrals in eq. ( |B1|) . Moreover, in this case, the BS vertex function Gm(po,p) depends 
upon the imaginary part of po, which allows one to use directly the numerical solutions 



obtained in the rotated system [2~1~,|2~0|. In particular, when (O) is defined at fixed value of 
p (which is just the case of eD processes investigated in this paper), the Wick rotation is 
no longer relevant in the computation of matrix elements. In order to establish a connection 
between the calculation of the matrix elements in the form (|Bl| ) and the matrix elements 
( p.30| ), we introduce the deuteron densities (see, also ref. IfTfl ) 



dp 1 



(2tt) (p 2 — in 2 + ie) 2 [j)\ — in 2 + ie) 
o E Tr {Gm(po, p)(pi + m)d(pi + m)G M (po, p)(p 2 + m)} , (B2) 



3 M 



where the integration contour is shown in Fig. 17. There are two kinds of singularities in eq. 



(|B2|) : i)when one of the particles is on mass shell, pio = w or P20 = w labeled "1" and "2" in 
Fig. |1^, and ii)when both particles are deeply virtual, pio = — w or p 2 o = —w , labeled "3" 
and "4", respectively. As previously mentioned, in eq. ( [B"2| ) the singularities are removed 
by performing the Wick rotation and by integrating along the imaginary axis of pq. Instead 
of rotating the contour, we close it in the upper semi-plane, and integrate eq. ( [B~2| ) in the 
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Minkowsky space. Thus there remain two singularities contributing to the full densities, 
a simple pole at p 2 o = w, and a pole of second order at p w = —w. The former, which 
corresponds to the spectator on the mass-shell, gives the dominant contribution to the full 
integral. It is exactly this contribution which enters our formulae in Section [Tl] (eq. ( |2.30| )), 
and in Section [TTI| (eqs. ( |2.39| )-( p.41| )). The approximation of the full matrix elements by 
the nucleon pole contribution is often used in the description of processes involving the 
deuteron Jl"3| . |3"2|j . Accordingly, let us define the nucleon pole contribution to the full density: 

(A\bs ( \ = _J 1 

\^/ P oi e yp) 2Md 2w m 2 d {m d - 2w) 2 

\Y. Jr {Gm(po, p)(pi + m)6(pi + m)G M (po, p)(p2 +m)\. (B3) 

It can be seen that the matrix elements ( (2.30|) is the pole part of ( pi|) . In so far as the pole 
contribution (eq. flB3|)) dominates the full density (eq. (|B^)), the relevant quantities (eqs. 
( p.39| )-( p.41| )) can be calculated by using the Wick rotation and by computing numerically 
the integral (|B2|) in the rotated system, using the numerical solutions of the BS equation 



obtained in ref. [^TJ. Fig. [18] shows the charge (70) and scalar (1) densities computed by eqs. 
( |B2| ) and ( |B3| ). It can be seen that up to |p| ~ 0.65GeV/c, both methods provide the same 
results (above 0.65GeV/c the pole contribution has been calculated using for the BS vertex 



function Gm, the analytical parameterization from ref. ||20|| . More details on the behaviour 



of various full and pole densities can be found in ref. [31 
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FIGURES 

FIG. 1. Diagrams corresponding to elastic eN scattering (a), and inclusive eD scattering in 
Impulse Approximation (b). 



FIG. 2. The diagram corresponding to the averaging of the operator O, eq. ( |2.20| ), between 
deuteron states. The crossed line corresponds to the nucleon on the mass-shell. 

FIG. 3. The Bethe-Salpeter matrix elements ( 1 }™f e (p) and (70 )^ Q f e (p) (solid lines) compared 
with their non relativistic limits, eq. ( |2.44| ), computed with the nucleon momentum distribution 
n £>(p) corresponding to the Bonn |^7j (dashed), Paris [^8) (dotted) and Reid j29| (dash-dotted) 
interactions. 



FIG. 4. The same as in Figj3| but for the vector density ( 73 )^,f e (p) 
corresponding to cos 9^ = 1. 



FIG. 5. The scaling function F (|q|, y), eq. vs |q| 2 for various values of y. For the sake 
of completeness, the value of the Bjorken variable XBj = Q 2 /2mu is also shown. 

FIG. 6. The dependence of the upper limit of integration |p| mai in eq. ( |2.35| ) upon |q|, for 
fixed values of y. 

FIG. 7. The asymptotic scaling function f BS (y) computed within the Bethe-Salpeter formalism 
(eq. (|3.5D evaluated at |q| — > 00). 



FIG. 8. Off-mass-shell contributions to / (y): the three curves correspond (from the top) to 
the three terms of eq. ( 2.41] ). 



FIG. 9. The "moving corrections" to F BS (\q\,y), eq. ( pT40| ) , for |q| = 3,10 and 18 GeV/c 
respectively. The solid line is the asymptotic scaling function f BS (y). 
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FIG. 10. The various contributions to the asymptotic scaling function f BS (y). Dot-dased 
line: the static part (|2.39| ); dotted line: the "moving" corrections ( [2.40 ). The solid line is the total 
scaling function (eq. ( |3.5| ) evaluated at |q| — > oo). 



FIG. 11. The BS asymptotic scaling function (full) compared with the non relativistic scaling 
function ( |3.2[) corresponding to the Reid (dashed), Paris (dotted) and Bonn (dot-dashed) interac- 



tions. 



FIG. 12. The Bethe- Salpeter S-wave | ^s(Po, |p|) I (solid line), with p = M D /2 - y/p 2 + m 2 , 
compared with the Gross .S-wave corresponding to the IIB solution |l3| (dot-dashed) and with the 
non relativistic deuteron S-wave obtained from the Paris (dotted) and Bonn (dashed) potentials. 



FIG. 13. The same as in Fig. 12 for the deuteron D-wave 



FIG. 14. The Bethe- Salpeter negative relative energy waves ^p 13 (po, |p|). 
FIG. 15. The Bethe-Salpeter covariant momentum distribution (eq. (|4.12| ) with 



Po = Md/2 — \/p 2 + m 2 ) (full) and the contribution from the negative relative energy states 
(|4l|) (dot-dashed). The non relativistic momentum distributions corresponding to the Bonn 
(dot-dot-dashed), Paris (dotted) and Reid (dashed) potentials are also shown. 

FIG. 16. The diagram defining the squared invariant matrix element for the deuteron in terms 



of the operator O, eq. (Al). The crossed line corresponds to the nucleon on the mass-shell. 



FIG. 17. The integration contour in eq. ( |B2[ ). The singularities labeled "1" and "2" correspond 
to the first and second particles on mass shell, pio = w and P20 = w, respectively, whereas 
the singularities labeled "3" and "4" correspond to both particles off mass shell, p±o = —w and 
P20 = ~w. 



•XA 



FIG. 18. The deuteron densities: a) (to)(|p|)- The solid line corresponds to the integration 
over po, and the dashed line to the nucleon pole contribution; the nucleon pole contribution to the 
scalar density (l)(|p|) is given by the dotted line, b) (l)(|p|). The solid line corresponds to the 
integration over po and the dashed line to the nucleon pole contribution; the dotted line represents 
the nucleon pole contribution to the scalar density calculated with (to)(|p|) (1 — P 2 /(^ m2 ))- 
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Fig. |2|. C. Ciofi degli Atti.... Relativistic Structure of the Deuteron. 



Fig. |j. C. Ciofi degli Atti.... Relativistic Structure of the Deuteron... 
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Fig. ||. C. Ciofi degli Atti.... Relativistic Structure of the Deuteron. 
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Fig. ||. C. Ciofi degli Atti.... Relativistic Structure of the Deuteron... 
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Fig. [7]. C. Ciofi degli Atti.... Relativistic Structure of the Deuteron... 
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C. Ciofi degli Atti.... Relativistic Structure of the Deuteron... 
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Fig. |9|. C. Ciofi degli Atti.... Relativistic Structure of the Deuteron... 
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10. C. Ciofi degli Atti.... Relativistic Structure of the Deuteron. 
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C. Ciofi degli Atti.... Relativistic Structure of the Deuteron. 
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Fig. 12. C. Ciofi degli Atti.... Relativistic Structure of the Deuteron. 
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C. Ciofi degli Atti.... Relativistic Structure of the Deuteron. 



47 



CO 




10- 4 



0,0 



0,5 1,0 
p, GeV/c 



14. C. Ciofi degli Atti.... Relativistic Structure of the Deuteron. 
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C. Ciofi degli Atti.... Relativistic Structure of the Deuteron. 
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C. Ciofi degli Atti.... Relativistic Structure of the Deuteron. 
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r ig. 17. C. Ciofi degli Atti.... Relativistic Structure of the Deuteron. 
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C. Ciofi degli Atti.... Relativistic Structure of the Deuteron. 



